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Abstract 
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1 Introduction 



Black holes, the most simple predictions of the general relativity, are one of the most 
enigmatic constructs in the present day physics. They draw physicists' attention as the 
paradigmatic objects to test possible quantum theories of gravity. Although still there is 
no convincing observational data in favor of conclusively proving the existence of black 
holes in the universe, there exist certainly sufficient evidences that make the study of such 
objects and the effects on their environment a matter of great importance to astrophysics. 
However, black holes are not objects of direct observing, so we must have to observe 
them indirectly through the effects they exert on their environment. They will greatly 
affect the surrounding plasma medium (which is highly magnetized) with their enormous 
gravitational fields, and hence, plasma physics in the vicinity of a black hole has become 
a subject of obvious interest in astrophysics. In the immediate vicinity of a black hole 
general relativity applies. It is therefore of interest to formulate plasma physics problems 
in the context of general relativity. 

The energy fiux carried out by black holes (or compact stars) produces a relatively 
large magnetic field. The study of stationary configurations and dynamic evolution of 
conducting fiuid in a magnetosphere of massive black holes demands the theory of general 
relativistic magnetohydrodynamics (GRMHD). The equations of GRMHD theory describe 
the aspects of interaction of relativistic gravity with plasma's magnetic field. The study 
of plasmas in the black hole environment is important because a successful study of waves 
will be of great value in aiding the observational identification of black hole candidates. 

An isolated black hole can have an electromagnetic field, if it is endowed with a net 
electric charge [H [21 [3l H] . Since a collapsed object can have a very strong effect on an 
electromagnetic field, it is of interest to determine this effect using GRMHD equations 
when a black hole is placed in an external electromagnetic field. 

A covariant formulation of the theory based on the fiuid equations in curved spacetime 
has so far proved unproductive because of the curvature of four-dimensional spacetime in 
the region surrounding a black hole. The 3+1 formulation of general relativity, developed 
by Thorne, et al O [6], [7], provides a method in which the electromagnetic equations 
and the plasma physics at least look somewhat similar to the usual formulations in fiat 
spacetime while taking accurate account of general relativistic effects such as curvature. 

In 3+1 formalism of Thorne, Price, and Macdonald (TPM) [8J, work connected with 
black holes has been facilitated by the replacement of the hole's event horizon with a 
membrane endowed with electric charge, electrical conductivity, and finite temperature 
and entropy. Mathematically the membrane paradigm is analogous to the standard, full 
general relativistic theory of black holes so far as physics outside the event horizon is 
concerned, and moreover, the formulation of all physics in this region turns out to be very 
much simpler than it would be using the standard covariant approach of general relativity. 

The pioneer work of 3+1 spacetime split, called ADM (Arnowitt, Deser, and Misner) 
formalism, was done in 1962 [9] to study the quantization of the gravitational field. Since 
then, their formulation has most been applied in studying numerical relativity [10]. TPM 
extended the ADM formalism to include electromagnetism and applied it to study elec- 
tromagnetic effects near the Kerr black hole. As a result, their work has opened up many 
possibilities for studying electromagnetic effects on plasmas in the black hole environment. 

In recent years there have been attempts to exploit the 3+1 formalism. Zhang [TT| [12] 
considered the case of perfect GRMHD waves in the vicinity of Kerr black hole and dis- 



2 



cussed the linearized waves for the cold (negligible particle pressure) plasma propagating 
in two- dimensions. Holcomb and Tajima [T3], Holcomb [13], and Dettmann et al. [15j 
investigated some properties of wave propagation in the Friedmann universe. Khanna 
[16] derived the GRMHD equations for two-fluid plasma in Kerr black hole. Anton et al. 
[17] investigated various test simulations and discussed magneto-rotational instability of 
accretion disks. Anile [18] worked on relativistic shocks/simple waves in magneto-fluids 
in cold relativistic plasma. Komissarov pTO] discussed the Blandford-Znajek monopole so- 
lution in black hole electrodynamics. Buzzi et al. [201 EI] described a general relativistic 
version of two-fluid plasma physics in TPM formulation and developed a linearized treat- 
ment of plasma waves in analogy with the special relativistic formulation of Sakai and 
Kawata [22] . They also investigated the one dimensional radial propagation of transverse 
and longitudinal waves near the Schwarzschild black hole. 

Recently, Sharif and Sheikh [23j investigated cold plasma in the vicinity of the Schwarzschild 
black hole horizon by using 3-1-1 formalism of the GRMHD equations and described the 
dispersion relation. 

In this paper we apply TPM formalism of the GRMHD equations to study the dy- 
namical magnetosphere of the Reissner-Nordstrom (RN) spacetime, which is the solution 
of the coupled Einstein-Maxwell equations and describes a spherically symmetric black 
hole endowed with electric or magnetic charge. The RN solution is important due to its 
intrinsic interest [2ll [25] specifically in the context of its recent applications in Hawking 
radiation [261 123 ISH [29] . If the horizon is sufficiently small, the magnetically charged RN 
solution develops a classical instability in the context of spontaneously broken gauge the- 
ories, which has significant implications for the evolution of a magnetically charged black 
hole [30]. It leads, in particular, to the possibility of evaporating a black hole completely, 
leaving in its place a nonsingular magnetic monopole. The magnetic monopole hypothesis 
was propounded by Dirac [311 132] relatively long ago. The ingenious suggestion by Dirac 
that magnetic monopole does exist in nature was neglected due to the failure to detect 
such a particle. However, in recent years the development of gauge theories has shed new 
light on it. Moreover, in asymptotically fiat space, the extremal RN black holes hold 
an important and controversial status in black hole physics. Usually, these black holes 
were supposed to be a limiting case of non-extremal black holes [33]. This conventional 
view was challenged in [3l] on the basis of the fact that the topology of the extremal and 
non-extremal black holes have qualitative differences. On the basis of these differences, it 
was argued that extremal RN black holes have zero entropy with no definite temperature 
in spite of having a nonzero horizon area [Ml [3S] . The extremal RN black holes are also 
important in the context of supergravity theories [361 [33, [38] • As shown in [39j , an exact 
solution of these black holes exists in (super) string theory. Thus aspects of the RN so- 
lution must be of interest in a broader contest. In view of this reason, our study of cold 
plasma in the environment close to the event horizon of the RN black hole is interesting. 
The result we have obtained reduces to that of the Schwarzschild black hole [23] when 
the charge term vanishes. 

This paper is arranged as follows. In Section 2, we summarize the GRMHD equations 
in the RN black hole magnetosphere in 3-1-1 formalism. We investigate the GRMHD 
equations for cold plasma in the case of rotating magnetized background in Section 3. In 
Section 4, we present our study of non-magnetized plasma in rotating background. We 
investigate the dispersion relations for the non-rotating background in Section 5. Finally, 
in Section 6 we present our remarks. We use natural units: G = c = 1. 
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2 3+1 Formulation of GRMHD Equations around 
RN Black Holes 



The 3+1 formulation of general relativity is based on the concept of selecting a preferred 
set of spacelike hypersurfaces which form the level surfaces of a congruence of timelike 
curves. The choice of a particular set of these hypersurfaces constitutes a time slicing of 
spacetime. The hypersurfaces considered here are of constant universal time t. 

In TPM notation, the metric of the Reissner- Nordstrom (RN) black hole with magnetic 
charges [iQl HI] is given by 

ds^ = g^^dx^'dx" = -fdt^ + ~ + ^^^^^^ + sin^ edif^), 
m = l- — + ^, Q' = Ql + Ql. (1) 



r 



The components denote spacetime coordinates and fi,i^ = 0,1,2,3. The vector po- 
tential of the dually charged RN black hole has non- vanishing components: At = Qe/^i 
A^p = —QmCOs9, with corresponding field strength: Fqi = Qe/r'^, -^23 = Qm/r'^, describing 
electric charge Qe and magnetic charge Qm- 

The metric ([T]) describes: (i) the Schwarzschild black hole for Q = 0; (ii) the mag- 
netically charged Reissner-Nordstrom black hole [3D] for Q^, = 0, < Qm < M; (iii) the 
generic Reissner-Nordstrom black hole for Qm = 0, < Qe < M; and (iv) the extremal 
Reissner-Nordstrom black hole for Q = M. 

According to the null hypersurface equation, the metric function /(r) can yield an 
event horizon at r = r_|_ and an Cauchy (inner) horizon at r = r_, where 

r± = M ± ^/^P'^Q^. (2) 

The curvature singularity at r = is hidden behind these horizons for Q < M, while 
for Q > M the metric ([T]) does not describe a black hole at all but rather a physically 
forbidden naked singularity, as no real value is found in ([2]). The mass parameter M of 
the RN solution is inseparably connected with the charge parameter Q [42j. This means 
that M = only when Q = 0. 

When Q = 0, r+ = 2M is the usual Schwarzschild horizon. In the extremal case of 
Q = M, the two values r± coincide and the horizon becomes degenerate. The extremal 
spacetime is just the flat space vacuum for Q = 0. An M > Q black hole will tend to 
Hawking radiate down to its extremal M = Q state. The Hawking temperature 



^ 27r(M+ v/M2-g2)2 

vanishes as the black hole attains its extremal limit. Thus the extremal black hole is 
typically quantum mechanically as well as classically a stable object. However, this is not 
true for extremal electrically charged black holes in our world [44]. The extremal black 
hole is quantum ground state of the charge Q superselection sector of the Hilbert space. 
This view is reinforced by the N = 2 supersymmetric version of the theory. 

The hypersurfaces of constant universal time t define an absolute three-dimensional 
space described by the metric ds'^ = Qijdx'^dx^ , where the indices i,j, refer to coordinates 
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in absolute space and range over 1, 2, 3. The fiducial observers (FIDOs), i.e. the observers 
remaining at rest with respect to this absolute space, measure their proper time r using 
clocks that they carry with them and make local measurements of physical quantities. 
Then all their measured quantities, such as velocities V and fields B and E are defined as 
FIDO locally measured quantities and all rates as measured by the FIDOs are measured 
using FIDO proper time. In making these measurements the FIDOs use a local Cartesian 
coordinate system that has basis vectors of unit length tangent to the coordinate lines: 

^'~^hr' rde' rsinO dip' ^ ^ 

The ratio of the rate of FIDO proper time to that of universal time is defined in terms 
of redshift factor: 

dr r 2M g2 

called the lapse function. It measures the amount of FIDO proper time elapsed during 
the passage of a unit amount of universal time. 

If a spacetimc viewpoint is considered rather than a 3+1 split of spacetime, the set of 
orthonormal vectors also includes the basis vector for the time coordinate given by 

The FIDO proper time r acts as a local laboratory time, where the FIDOs have the role 
of "local laboratories". It does not provide a slicing of spacetime as it is not a global 
coordinate. The RN time coordinate t is the logical choice to satisfy this role and in fact 
slices spacetime in the way that the FIDOs would do physically. All subsequent equations 
are therefore expressed in terms of the universal time coordinate t rather than the FIDO 
proper time r. 

The lapse function a acts as a gravitational potential and governs the ticking rates 
of clocks and redshifts as well. It also measures the gravitational acceleration felt by a 
FIDO: 

a = Vlua = i('^-^)e,. (7) 

The rate of change of any scalar physical quantity or any three-dimensional vector or 
tensor, as measured by a FIDO, is defined by the convective derivative 

" ^i^.V.vV (8, 



Dt \adt 

where V is the velocity of a fluid as measured locally by a FIDO. 

The Rindlcr coordinate system, in which space is locally Cartesian, provides a good 
approximation to the RN metric near the event horizon. However, the essential features 
of the horizon and the 3+1 split are retained without the complication of explicitly curved 
spatial three-geometries. The RN metric is approximated in Rindler coordinates by 

ds'^ = -a'^dt'^ + dx'^ + dy^ + dz'^, (9) 

where 



X 



= - 1) , y = r+ip, z = 2r+a. (10) 
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The standard lapse function is again denoted by a and simplifies in Rindler coordinates to 
2;/(2r_|_) , where r_|_ is the location of the event horizon of the RN black hole. This function 
vanishes at the horizon which we can place at z = and it increases monotonically as z 
increases from to oo. 

Maxwell's equations in 3+1 formalism take the following form: 

V-B = 0, (11) 
V-E = 47rp„ (12) 

^ = -Vx(aE), (13) 
dE 

— = V X (aB) - 47raj, (14) 

where pe and j are electric charge and current density, respectively. For the perfect MHD 
(i.e., MHD with perfectly conducting) assumption there exists no electric field in the 
fluid's rest frame, i.e., E + V x B = 0. Under this condition the equation for the evolution 
of magnetic field (|T3|) becomes 

_ = V X («V X B) (15) 
The conservation of mass and momentum equations are written, respectively, as follows 

^ + {(aV).V}(p„/i) 
(9V 

+p„/i7'V ■ — + poP^^V ■ (aV ■ V)V + pop{V ■ (aV)} = 0, (16) 



1 9 „\ fB\ 1 

1 



= -pWa, - p,i + -^(V X B)jV . (V X B) - ^^(aB) 

+ ^(aB.^^,B^ - J-[B X {V X (V X (aV x B))}]„ (17) 

where a subscript i on a vector quantity refers to the i-component of that vector. The 
p = {p + p) / Po is the specific enthalpy of the fluid, where p is the total density of mass- 
energy and p, the pressure as seen in the fluid's rest frame. The po is the fluid's rest-mass 
density and 7 = (1 — V^)~^/^ is the fluid's Lorentz factor as seen by the FIDOs. Equations 
([I5])-([I7]) are the perfect GRMHD equations for the RN black hole. 

We assume for simplicity that the plasma has vanishing thermal pressure and vanishing 
thermal energy density; that is, we restrict our investigation to a "cold plasma," for which 
we consider that the total density of mass-energy p remains the same as the rest-mass 
density po'- 

p -\~ p 

p = 0, p = Po, and p = = 1. (18) 

Po 
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Now using (fT8|) in (fT5!) - (fT7j) we obtain the perfect GRMHD equations for cold plasma 
near to the event horizon of RN black hole. 

We characterize the perturbed flow in the magnetosphere by its velocity V and mag- 
netic field B as measured by the FIDOs, and the fluid's density p. The first order per- 
turbations in these quantities are denoted by 5V, 5B and 5 p. Accordingly, the perturbed 
variables take the following form: 

B = B° + (5B, V = V° + 5V, p = p° + 5p, (19) 

where B°, V° and p° are unperturbed quantities. The waves can propagate in z-direction 
due to gravitation with respect to time t and thus perturbed quantities must depend on 
z and t. 

3 Cold Plasma in Rotating Magnetized Background 

We use the linear perturbation and Fourier analyze techniques to reduce GRMHD equa- 
tions to ordinary differential equations. The magnetosphere has the perturbed flow along 
x-z plane in this background. The FIDO-measured fluid four-velocity can be described 
in this plane by 

Y = V{z)e^ + u{z)e^, (20) 
while the Lorentz factor 7 takes the form 

^ ~ V1-U2^W ^^^^ 

The rotating magnetic field can be expressed in the x-z plane as 

B = 5[A(;2)e, + e,]. (22) 
The variables A , u and V are related by 

V = ^ + \u, (23) 
a 

where Vp is an integration constant. 

We use the following perturbation quantities 

p = 5p/p = p{t,z), V = 5V = Vx{t,z)e,^ + V;,{t,z)e^, b = ^ = h^{t,z)ej, + h^{t,z)e:„ 

(24) 

with 

p(t, z) = cie-*('^*-'=^), v,{t, z) = cae-^^-*-'^^), v^it, z) = c^e-'^^'-^^^ 

Kit, z) = C4e-*('^*-^'^), 6,(t, z) = cse-^^^*-'^"), (25) 

where c^, s = 1, ■ ■ ■ , 5, are arbitrary constants. 
Using dMD in (US}- (HZ]), we obtain 

= V X (av X B) + V X (aV X (5B), (26) 
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V ■ (SB) = 0, 



aot J \aot J P 

= -2p-f\Y ■ v)(V • V)ln7 - p-f^(V ■ VV) ■ v + p(v -Vlnu), 



2 , 



1 



1 



47r/"^^ ' ^' 47r^ M a dt ' A-n 



1 1 
+ — 



B X <^ V X 



1 d{5B) 
a dt 



1 



= -7'{(5p + 2p72(V • v)}a, + - {aBj),,}SB^ 

Ana 

-pi\viV' + v'Vi)Vk,jV^ - -i''{5pV' + 2p7'(V ■ w)V^ + pi;^'}\^.j 
-7^yi{5p\/^' + 4p72(V ■ w)V^ + pv^}Vj^kV^. 

The component form of fl26|) - fl29|) can be written as follows : 
1 dh 



a dt 
a dt 



+ ubx 



+ ubz 



V In a{vx — Xvz + Vb^ — ub^) 
+ - Af2,^ - X'vz + V'bz - u'b^), 

0, 
0, 



1 dp 



a at \adt 



1 (9u^ 



+ M^^a;,^ + 7 M 



1^ 

a dt 



+ (1 + 7V)t;, 



-7'm{(1 + 2YV')V' + 2YuVu}v, 



u 



P72(1 + 7V2) + — ^ — P7^MV^- 



+<'p7"(l+7V^ 



A-n ] a dt 



, I M^^x,^ + P7 «^ + , 

47r \ 47r 



An J a dt 



B 



uvz,z - -;— (1 - u % 
An 



B' 



Ana 



b^{a'{l - u^) - auu'} + ppYu{{l + YV^)V' + Y^u'V} 



+ 



P7^m{(1 + 47V^)W + 4(1 + YV')VV'} + 



B^ua' 

Ana 



+ [p7^ [{(1 + 2^\?){1 + 27V2) - 7Vn^' 



B\ 



+27^(1 + 27V)WVl + =^(aXy 

Ana 



Vz = 0, 
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P7 (1 + 7 M ) + — — ) — —+ uV - 



+ P7^(1 + 7V) 



47r j a dt 
A^B^ 



An / a dt 



An 



uvz,z + [ p'y uV + 



AB^ 

47r 



AB^ g2 

+ — — (1 — u^)bx z + - — {{oiXy + a'X — u\{ua' + u'a)}bx 
An ' Ana 

+p-f^[az + u{{l + -f\^)u' + -f\VV'}] 



+ 



XB'^ua' 
Ana 



P7^{uV(l + 47V) + 2V{az + uu{l + 27%'))} + 



+ [P72{m'(1 + 7V)(1 + 47 V) 
+2u72{(l + 2-f^u^)VV' + a,}} - :i^(aA)>, = 0. 



Ana 

From the Fourier analyzed of fl30|) - fl35|) with fl25|) . we obtain 

C3(a' + ika) — C2{(a;A)' + ikaX} + 05(0;^^)' — C4{{au)' — iu + ikau] 

( iuj \ 

C5 h iku = 0, 

\ a J 

c^ik = 0, 



+ci ( — — + iku ] + C2 

a 



7^M + iA;(l + 7^^) 



a 



u 



(1 - 2^^v?){\ + ^^v?)- + 27^uVl^' 



u 



+C37 



+ ifcu M/ + u{(\ + l^'V'W' + 27%mV} 

a 



CiP7\{(l + 7V2)\/' + fm^V} - ( (1 - + -(1 - _ 

477 a 



+C2 



— ( p^'^uV — ] V iku ip^^vV + 



477 / a 

+P72{(1 + 1^''u^){\ + 27V) - 7V}\/' 

V?u 

+2p7^(l + 272m2)mmV + -— (Aa)' 



477 



477a 



+ C3 



P7'(1+7V2) + |^ 



B 

+zA;m<{p72(1 + 7V2)- — 



+P7^m{(1 + 47V2)W + 4(1 + -fV^)VY} 



B^ua'' 
Ana 



0, 



+C2 



P7^(l +7V) + 



An 



+ <'P7'(1 + 7V)- 



+ <; P7'K(1 + 7'«')(1 + 47'^^') + 2V((1 + 27V)\/F' + a,)} - ^^(aA)' 



C3 



— ( p'j'^uV — 



XB'\ iuj 



An I a 



AB^ 

47r 



+ <^ p7'{w'^'(l + 47V2) + 2V{a, + uu'{l + 27^2))} + 



A(l - u^)tk + A(l - u^)- - Xuu' + 



AB^wa' 

Ana 



+— C4 
47r 



(41) 



Equation (1381) implies that C5 is zero which gives that bz is zero. The determinant of 
the coefficients of ci, C2, C3 and C4 in (l36ll and (l39l) - (HTl) equated to zero gives a complex 
relation in k, called the dispersion relation, which is of the form A[z,uj)k^ + B[z,uj)k^ + 



Numerical Solution Modes 

We investigate the different types of modes of waves when B > and the wave number 



is in arbitrary direction to B. We use the lapse function a 



2r+ 



< < 1- We consider the black hole mass M ~ IMq, Q^/M^ 

47r 



2. From the mass conservation law in three-dimensions we get u 



= 0.7, p = 1 and 
For 



simplicity, we also assume that u = V. From fl23l) we get A = 1 — 



V2H 



by taking Vp 



which shows that the magnetic field diverges close to the horizon. 

Using these values in the dispersion relation we get values for k, from which we have 



calculated the phase velocity Vp = f and group velocity Vg 



[n+LU 



dn\ — l 



where n(= 1/fp) 



is the refractive index computed as the ratio of the speed of light in a vacuum to the speed 



of light through the material and ^ determines whether the dispersion is normal or not. 

The real part of the dispersion relation gives four real values for k , out of which two 
are real and interesting. The other two values turn out to be imaginary in the whole 
region. One real value for k is obtained out of three values found from imaginary part of 
dispersion relation. This is shown by Fig. 3. The two dispersion relations obtained from 
the real part are shown in the Fig. 1 and Fig. 2. 

We see from Fig. 1 that the waves gain energy with the increase in angular frequency 
but lose when we move away from the horizon and hence damping arises. But in the 
vicinity of the horizon the wave number is very large due to strong gravitational field, 
so no wave exists there. The phase and group velocities increase as we depart from the 
horizon. The dispersion is not normal there because of refractive index n < 1 though 
^ > 0. 

did 

The Fig. 2 indicates that the wave number is proportional to angular frequency and 
inversely proportional to z. There exists no wave at the horizon since wave number is 
infinite there. The phase and group velocities show the same behavior. We observe the 
damping modes as we depart from the event horizon and the growing modes as z decreases. 
Since n <1 and ^ < 0, the wave does not disperse normally there. 
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Figure 3 shows that the wave number is very large at the event horizon {z = 0) but it 
decreases as we depart from the horizon and also it takes negative value for some values 
of z and u. The phase and group velocities are also negative and the group velocity is 
greater than the phase velocity. The value of refractive index n < 1 and ^ < for some 
values of z and u. Therefore the region is of anomalous dispersion and possesses the 
properties of metamaterials. 



4 Cold Plasma in Rotating Non-magnetized Back- 
ground 



In non-magnetized background B = i.e. when fo^; 
(l36ll - (HTl) reduce to the following form: 



0, the GRMHD equations 



Cl 



a 



+ iku + C2 



+ (1 + 7V)zA; - (1 - 2^\^){1 + ^\^)- + 2^\^VV' 



a 



u 



V + ikuY + 7'm{(1 + 27V)\/' + 27'uuV} 

OL 



0, 



(42) 



h iku 7 uV 

a. 



Ci72u{(l + 7V2)l^' + 7\uV} + C27^ 
+ {(1 + 27V)(1 + 27^2) - 7V2}1^' + 27^(1 + 2^''u^)m^V\ 



+C3 



a 



+ iku 7^(1 + 7V 



+7^M {(1 + ^'fV'^)uu + WV\\ + 7V^)}] = 0, 



(43) 



Ci7^{a2 + (1 + ^'^u^)uu + 7WV"} + C2[7^(l + 7V) (^-^ + i^^^ 
^f{u!{\ + 7'u2)(l + 472^2) + 2u72(a, + (1 + 2^''u^)VV')}\ 
+C37V {-^ + ik^ u + u^{\ + 47V2) + 2{a, + W(l + 272^2)} 



0. (44) 



where the FIDO-measured fluid four-velocity V, Lorentz factor 7 are given by fl20l) and 
fl2T|) . respectively. The determinant of the coefficients Ci, C2 and C3 in p2|) - fH^ yields a 
complex dispersion relation of the form A{z,ijj)k^ -|- B{z,uj)k'^ + C{z,uj)k + D{z,u) = 0. 

To analyze the numerical solution mode we consider Q'^/M'^ = 0.5 and assume that 
V = u and p =constant. The mass conservation law gives that u = r-^-^ . From the real 

' " Vz +2 

part of the dispersion relation we get only two real values for k which are shown in Fig. 4 
and Fig. 5. The imaginary part gives only one real root and two complex conjugate roots. 
Figure 6 shows the real value for k obtained from the imaginary part of the dispersion 
equation. 

We see from Fig. 4 that the wave number is very large close to the event horizon and 
the waves lose energy as we move away from the horizon of RN black hole. This shows 
that the increase in u increases k and the waves are in growing mode as z decreases. The 
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phase and group velocities are of the same pattern and increase as z increases. Since 
n < 1 and ^ < 0, the region is not of normal dispersion. 

Figure 5 shows that the real waves exist in the region but no wave in the vicinity of 
the event horizon because of the existence of strong gravitational field there. The wave 
number increases as angular frequency increases but decreases as z increases. The wave 
number decreases as we go away from the event horizon and hence damping takes place. 
For this region refractive index n < \ and ^ > which implies that the region is of 
anomalous dispersion. 

In Fig. 6, infinite wave number occurs at the horizon, so no wave exists there. The 
phase and group velocities admit the same pattern. These velocities increase when we 
depart from horizon. The region has n < 1 and ^ < for some values of z and uj. Hence 
it is not the case of normal dispersion. 

5 Cold Plasma in Non-rotating Background 

The magnetosphere has the perturbed flow only along z-axis in this background. Hence, 
FIDO measured fluid four-velocity V is described by 

V = u{z)&^. (45) 

The Lorentz factor takes the form 

V 1 — M 

while the FIDO measured magnetic field becomes 

B = 5(^)e,. (47) 
Then the equations (!36|) - PT|) reduce to the following form : 



a 



C5 = (4^ 



ifccs = (49) 

(Z^^,^k^ +c2|(1 + 7V)2A;-(1-27V)(1 + 7V)---7'm1 =0 (50) 
\ Ci ) \ u a ) 

cxl^{az + uv!{\ + 7V)} + C2[7^(l + 7^M^) + i'^u^j 

+72{m'(1 + 7V)(1 + 47V) + 2u72aJ] = 0. (51) 

From fH5]) or fH^ . C5 = 0; hence, there does not exist any perturbation in magnetic field 
of the fluid. We also get the same two equations fl5U]) and fl^ for the non-magnetized 
plasma. The complex dispersion relation, that follows from (!50l) and (15T!) . is 

mV(1 + nV)^^ - {m^7^(-l + mV) + ^^(1 + nV)(3iu^uV + — )}fc 



a 

7 { — (a^ + MM + U u 7 )(1 — M 7 — 2^ 7 ) 



u 

2 

-2{M' + M(a, + 4W)72 + 3MV7n--(l + MV)^} =0. (52) 
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We investigate the longitudinal waves propagating parallel to the magnetic field B. 
We analyze the numerical modes for Q^/M^ = 0.5. From the mass conservation law we 
get u = ^^-^-j-j- . The real part of fl52l) gives us two values for k in terms of z and u, shown 
in Fig. 7 and Fig. 8, while the imaginary part provides one real value for k, shown in 
Fig. 9. 

Figure 7 shows that the wave number k decreases as z increases i.e. the waves are 
growing energy with increase in u and decrease in z but damping occurs when z raises. So 
waves drop energy when we depart from event horizon. At the event horizon of the hole, 
we see that the wave number becomes infinite which means that the waves disappear due 
to the effect of immense gravity. The group and phase velocities are same. The refractive 
index n > 1 and ^ > in this case. Hence the dispersion is normal. 

We observe from Fig. 8 that the wave number is infinite at 2; = and hence no wave 
exists there. The wave number decreases as we depart from the event horizon. The waves 
show damping modes for increasing z. The increase in u increases k. The phase and group 
velocities have the same behavior. Since ^ < 0, the region is not of normal dispersion. 

Figure 9 shows that the wave number becomes very large and hence there exists no 
wave very close to event horizon. The waves are gaining energy with the increase in u but 
losing with the increase in the distance from the event horizon. The surfaces of group and 
phase velocities are showing the same nature. Since n < 1 and ^ < for most region, so 
the region is not of normal dispersion. 



6 Concluding Remarks 

Our main concern has been the investigation of the wave properties of cold plasmas in 
the RN black-hole's magnetosphere. We have exploited the 3+1 formalism of general 
relativity, developed by Thorne et al [H [6l [71 [8] , and derived the GRMHD equations in 
component form by linear perturbation. These equations are then Fourier analyzed to ob- 
tain dispersion relations. We have considered both non-rotating and rotating backgrounds 
(either non- magnetized or magnetized). The properties of plasma waves are analyzed on 
the basis of the quantities: wave number, phase and group velocities, and refractive index, 
derived from the dispersion relations. We solve the dispersion relation numerically and 
present the results graphically. 

Our investigation shows that the wave number becomes infinite at the event horizon 
and consequently, no wave is present there due to immense gravitational field. This 
indicates that no signal can pass the event horizon or near to it. But when we depart 
from horizon, the waves lose energy. Wave number is directly proportional to angular 
frequency and inversely proportional to z. Therefore the waves are in damping mode as 
we go away from the horizon and in growing mode as we approach the horizon. 

For a non-rotating background, the magnetospheric fiuid disperses normally. Figure 

7 shows the normal dispersion, while Fig. 8 and Fig. 9 show dispersions which are not 
normal. In the case of a rotating non-magnetized background, shown in Fig. 4, Fig. 5, 
and Fig. 6, we have found cases which are not normally dispersive. 

For a rotating magnetized background, we have found a region, shown in Fig. 3, where 
wave number, phase velocity, group velocity, refractive index, and ^ are all negative. This 
implies that the region has all the characteristics of left-handed metamaterials. The group 
velocity is greater than the phase velocity and 4^ < 0, which mean that the region is of 
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anomalous dispersion. The dispersion relations provide that the region under discussion 
is of non-normal dispersion. Hence real signals cannot pass through this region. 

Our study shows that the wave number becomes infinite at the event horizon and hence 
no wave exists there, which supports the well-known point of view that no information 
can be extracted from a black hole. Mackay et al. [15] found that rotation of a black hole 
is required for the existence of negative phase velocity propagation and the waves of less 
angular velocity are evanescent. It is interesting to mention that our analysis shows that 
negative phase velocity propagates in the rotating background whether the black hole is 
rotating or non-rotating. 

The MHD waves are non- dispersive in the cold plasma, but our analysis, shown graph- 
ically, predicts that they are dispersive. This happens because of the formalism used and 
the equations which are different from the usual MHD equations. Since the 3-1-1 split 
of general relativity is used in investigating waves propagating in a plasma influenced by 
the gravitational fleld, internal gravity waves which interrupt the MHD waves imply the 
cases of dispersion in each of the hyperfsurface. The flgures of the waves are given in a 
particular hypersurface of constant time t but not in the whole RN background. Hence 
our flnding is justifled locally, not globally. 

The result, obtained in this paper, reduces to the case of the Schwarzschild black 
for Q = 0, as was obtained in [23] ■ As described in the introduction, the Reissner- 
Nordstrom solution describes several special charged solutions in special cases, such as, 
(a) the magnetically charged solution for Qe = 0, < Qm < M, (b) the generic solution 
for Qm = 0, < Qe < M, and (c) the extremal solution for Q = M. Thus the result of 
this study can be specialized for these interesting solutions by suitably choosing the black 
hole parameters. The extreme RN black hole is distinguished by its coldness (vanishing 
Hawking temperature) and its supersymmetry. It occupies a special position among the 
black-hole solutions to the Einstein or Einstein-Maxwell equations because of its complete 
stability with respect to both classical and quantum processes permitting its interpretation 
as a soliton [1^ Wl\ . The extremal spacetime is also special in admitting supersymmetry 
in the context of = 2 supergravity [13 SSI SHI [SHI El [S2]- Thus, aspects of the 
Reissner- Nordstrom solution might be of interest in a broader context. 

In view of the above reasons, our study of the dispersion relation for the cold plasma 
near the event horizon of the RN black hole is well motivating. The results of this paper 
can be extended to the RN spacetime generalized with a cosmological parameter. This 
type of extension may be interesting from the point of view of an inflationary scenario of 
early universe. 
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Surface plot of wave number 




Surface plot of phase velocity Surface plot of group velocity 




Figure 1: The region is of non- normal dispersion, phase and group velocities are the same 
and increasing with increase oi z, > but n < 1. 
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Surface plot of wave number 




Figure 2: The region is of non- normal dispersion, phase and group velocities are the same, 
n < 1 and ^ < 0. 
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Surface plot of wave number 




Figure 3: The region is of anomalous dispersion, Vg > Vp, and k,Vp,Vg are negative, 
n < 1, 4^ < 0. The medium has the properties of metamaterials. 
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Surface plot 



of wave number 
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Surface plot 



of wave number 




Figure 5: The region is not of normal dispersion, phase and group velocities are the same, 
n < 1 and ^ > 0. 



21 



Surface plot of wave number 
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Surface plot of wave number 
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Surface plot of wave number 
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Surface plot of wave number 




Figure 9: The region is not of normal dispersion, phase and group velocities are the same, 
n < 1 and ^ < in most region. 
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